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Application  of  Mixed  Hj/EL  Optimization 
Abstract 

A  nonconservative  optimization  technique  is  utilized  to 
examine  the  problem  of  minimizing  the  2  norm  of  one  transfer 
function  subject  to  an  oo-norm  bound  on  another  transfer 
function.  For  single  exogenous  input /single  exogenous 
output  systems  this  nonconservative  mixed  solution  is  shown 
to  be  up  to  5%  better  than  a  previous  technique  in  a 
specific  example.  For  single  exogenous  input/two  exogenous 
output  systems  the  mixed  solution  is  shown  to  be  up  to  27% 
better  than  a  previous  technique  for  a  specific  example. 
Finally,  a  two  exogenous  input /two  exogenous  output 
practical  system  is  examined  to  demonstrate  the  utility  of 
this  mixed  optimization  technique. 
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Application  of  Mixed  Hj/H.  Optimization 


Chapter  I.  Introduction 


1 . 1 '  Background 

"The  designer  of  a  Single-Input-Single-Output  (SISO)  or 
Multiple-Input-Multiple-Output  (MIMO)  system  must  contend 
with  two  major  factors  in  the  design  process.  First,  the 
design  must  minimize  the  energy  of  the  output  errors  in  the 
face  of  input  noises.  Second,  the  design  must  minimize  the 
output  errors  given  a  norr-noise  bounded  energy  input  that 
would  affect  robustness.  Noise  sensitivity  will  be  called 
performance  in  this  work,  and  using  the  LQG  technigues  in 
modern  control  theory  is  easily  handled.  LQG  by  itself 
however  usually  exacerbates  robustness  problems.  LQG/LTR 
addresses  robustness  in  an  unstructured  manner,  but  an 
acknowledged  shortcoming  in  modern  control  design  has  been 
the  lack  of  visibility  into  performance /robustness  trade¬ 
offs  . 

Recent  efforts  have  focused  on  characterizing  this 
problem  with  the  2 -norm  as  a  measure  of  performance  and  the 
00-norm  as  a  measure  of  robustness.  These  are  usually 
referred  to  as  Hj/H*  Optimization  problems  and  they  can  be 
formulated  in  a  variety  of  different  ways.  The  formulation 
for  this  thesis  can  be  stated  as  maximizing  performance  at  a 
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given  robustness  level.  Numerous  approaches  to  this  problem 
currently  exist,  but  only  one  is  non-conservative.  The  non¬ 
conservative  approach  was  first  forwarded  in  1991  by  Ridgely 
[1]  and  will  be  referred  to  as  the  Mixed  Solution  throughout 
this  thesis. 

The  Mixed  Solution  has  only  recently  become  available 
and  allows  Hj/H,  trade-offs  to  be  made  for  a  wide  class  of 
systems.  The  robustness  design  parameter  ;an  be  used  to 
handle  a  variety  of  problems  such  as  unmodeled  dynamics, 
sensor  bias  and  disturbances,  excessive  gain  scheduling,  and 
input  disturbances. 

1.2  Purpose 

The  purpose  of  this  thesis  is  to  demonstrate  the  use  of 
the  Mixed  Solution  and  compare  results  to  other  techniques. 
First,  the  Mixed  Solution  will  be  compared  to  the  Central 
controller  for  a  system  with  a  single  exogenous  input  and  a 
single  exogenous  output  (  I 10).  Then,  the  Mixed  Solution 
will  be  compared  to  the  technique  forwarded  by  Bernstein  and 
Haddad  [5]  for  a  system  with  one  exogenous  input  and  two 
exogenous  outputs  (1120).  Finally,  a  physically  motivated 
system  with  two  exogenous  inputs  and  two  exogenous  outputs 
(2120)  will  be  examined  using  the  Mixed  Solution  to 
demonstrate  the  practical  aspects  of  this  technique  (no 
other  two-in  two-out  techniques  exist  for  comparison). 
Various  plants  will  be  used. 
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1.3  Overview 


In  Chapter  2,  various  aspects  of  control  theory  and  the 
role  of  Hj  and  theory  are  discussed.  In  Chapter  3,  the 
procedure  and  analysis  techniques  required  to  obtain  a  Mixed 
Solution  are  discussed.  In  Chapter  4,  problem  synthesis, 
description,  and  results  for  the  1110  system  are  discussed. 
Likewise  in  Chapter  5,  the  1120  system  is  covered  and  in 
Chapter  6  the  2120  system.  Chapter  7  contains  a  summary  and 
conclusions  followed  by  recommendations  for  further  study. 

The  software  package  PRO-MATLAB  was  used  for  this 
thesis.  FORTRAN  77  (UNIX)  was  also  utilized. 
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Chapter  II.  Mixed  Hj/H,  Control  Theory 


2.1  Singular  Values.  H-,  Norm,  and  Norm 

The  basic  analysis  tool  in  multivariable  control  design 
is  the  singular  value  versus  frequency  plot.  The  singular 
values  of  a  matrix  A  are  related  to  the  eigenvalues  of  A  as 
follows 


o,(A)  -  Ik^  (A‘A) 


(2.1) 


with  tJi(A)  =  ith  singular  value 

Xi  =  ith  eigenvalue 

A*  =  complex  conjugate  transpose  of  A 
A  e  C"" 

The  largest  singular  value  is  denoted  o  and  the  smallest  is 
denoted  o.  These  two  scalars  are  a  measure  of  the  size 
(gain)  of  a  transfer  function  matrix  at  a  given  frequency. 
Thus,  o  and  o  are  frequency  dependent.  In  SISO  systems,  o  = 
o  and  a  plot  of  a  verses  frequency  is  just  the  Bode 
magnituie  plot.  The  2-norm  is  related  to  the  area  under  the 
o  vs  CO  (frequency)  curve  and  is  given  by 

||G(  jw)  j  .  Lr  IG*(  j(i))  G(jbj)  ]  dtJ  (2*2) 

By  Parseval's  Theorem  it  can  be  seen  that  this  is  the 
transfer  function's  total  energy  output.  This  nicely 
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characterizes  the  energy  of  a  system  in  the  presence  of 
white  noise. 

The  oc-norm  of  a  transfer  function  is  given  by 

[G{ju))l,^supa  iGijw)]  (2.3) 

(i) 

and  is  a  measure  of  the  maximum  gain  of  the  transfer 
function.  Thus,  a  singular  value  plot  with  a  large  spike, 
and  hence  large  oc-norm,  would  not  be  considered  robust  at 
the  spike's  frequency  because  any  excitation  at  the 
frequency  of  the  spike  causes  a  potentially  disastrous  gain. 
Another  way  of  looking  at  the  oc-norm  this  is  through  the 
small  gain  theorem.  Suppose  we  have 


Figure  2.1  Small  Gain  Block  Diagram 

with  G  and  A  E  RH,.  Then  this  system  is  closed  loop  stable 
if 


lAIL-llGlL  <  1 


(2.4) 
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If  IGI,  <  Y  then 


lAI.  <  ^ 


(2.5) 


is  sufficient  to  guarantee  stability  of  the  system.  Thus, 
the  smaller  y  (JG|,)  is,  the  larger  |A|,  can  be  and  still 
have  guaranteed  system  stability.  If  we  think  of  A  as  a 
system  disturbance  or  uncertainty  and  G  as  a  nominal  closed- 
loop  system,  then  keeping  y  low  will  maximize  the  system's 
disturbance  tolerance  or  robustness. 

2 . 2  Motivation  and  Characterization 

Any  physical  system  will  have  noise  and  uncertainty  in 
it.  Thus,  it  makes  sense  to  want  to  minimize  the  system's 
sensitivity  to  noise  and  maximize  its  robustness.  In  this 
thesis,  this  problem,  called  mixed  Hj/H,  optimization,  is 
formulated  as; 

infimize  I T^^,#  j  subject  to  the  constraint 

K  stabilizing  that  iT^^jl,  s  y 

With  a  system  as  shown  in  Figure  2.2,  T^^,  indicates  the  z  to 
w  transfer  function  while  T^  the  e  to  d  transfer  function. 
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d 


c 


Figure  2.2  Standard  Problem  Block  Diagram 

The  w  input  will  normally  be  white  noise  while  the  d 
input  will  be  bounded  energy  disturbances.  The  following 
definitions  are  helpful: 

Yo  =  inf  iTedll, 

K  stabilizing 

=  inf  |T,„|lj 
K  stabilizing 

Kjopt  “  K  required  to  achieve  oc,, 

Yz  =  lT«d!.  when  K  = 

For  clarity  in  the  following  discussion,  consider  T„ 
T^.  First,  it  can  be  seen  that  minimizing  |  T^^l  ^  will 
produce  a  system  of  finite  (and  probably  low)  bandwidth  in 
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order  to  minimize  the  area  under  the  plot.  However, 
minimizing  iTe^U  typically  produces  a  closed-loop  system  of 
infinite  bandwidth,  as  the  low  frequency  energy  is  forced  to 
higher  and  higher  frequencies.  Thus,  mixed  Hj/H, 
optimization  has  competing  objectives. 

When  BTjwIj  is  at  the  resulting  controller  (K^opt)  is 
unique.  Thus,  Kj^p^  will  produce  a  finite  iT^dl.  {Y2)  which  is 
as  large  as  will  ever  get  in  the  mixed  H;/E.  problem. 

As  |Tedj«  approaches  y^,  |Tj„|lj  goes  to  infinity  since  the 
closed-loop  bandwidth  (and  thus  area  under  the  curve) 
typically  becomes  infinite.  A  typical  trade-off  would  look 
like  that  in  Figure  2.3. 

No  matter  how  large  y  gets,  iTg^l,  will  never  exceed  y^. 
In  fact,  it  really  makes  no  sense  to  talk  about  mixed  Hj/H. 
solutions  with  y>Y2  since  there  is  no  trade-off,  and  the 
solution  is  just  trivially  Kj^p^.  The  above  plot  has  been 
shown  [1]  to  be  monotonically  decreasing  for  the  Mixed 
Solution . 

Of  key  importance  is  that,  for  Y>Yof  the  controller 
that  gives  (T^U  «  Y  is  not  unique.  An  infinite  number  of 
stabilizing  controllers  can  achieve  iT^dl,  s  y-  Thus,  for 
Yo<Y<Y2/  the  fcimily  of  stabilizing  controllers  will  contain 
a  controller  such  that  iT^wlj  is  a  minimum  given  iT^dl-  ^  Y* 

In  (1:120)  it  is  shown  that  this  minimum  is  achieved  with  a 
controller  such  that  =  Y*  The  curve  given  by  a  plot 

such  as  in  Figure  2.3  is,  for  the  Mixed  Solution,  a  true 
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controller  or  a  full  state  feedback  controller.  The 
exogenous  inputs  in  the  d  vector  represent  conunands, 
disturbances,  etc.  into  the  plant  P.  The  regulated  output 
in  the  e  vector  is  the  system  response  that  will  be 
associated  with  the  oo-norm.  The  exogenous  inputs  in  the  w 
vector  represent  white  noises  coming  into  the  plant  P.  The 
regulated  output  in  the  z  vector  is  the  system  response  that 
will  be  associated  with  the  2-norm.  The  output  y  is  fed  to 
the  controller  K  which  generates  the  u  input  into  P.  The 
plant  can  be  represented  in  partitioned  transfer  function 
form  as 


(2.6) 


which  yields: 


e  =  P^d  +  P^^w  *  P^^u 

z  =  ^  P^^w  - 

K  =  P^rfd  -  P^^w  *  Py^u 


Using  a  lower  linear  fractional  transformation,  the  transfer 
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functions  T„  and  are: 


"  P.ui<  ^^-PyJ(^  ’  Py. 
Ped  =  Pod  ^  P.uf<  U-P,,K]  ’  Pyrf 


Note  that  the  characteristic  equation  and  hence  system 
stability  is  determined  by  the  [I-Pj^K]"^  term.  There  are,  of 
course,  "cross  transfer  functions"  and  but  neither 

of  them  is  directly  addressed  in  the  mixed  problem.  They 
really  don't  need  to  be  included,  since  any  concern  about 
these  cross  transfer  functions  should  be  taken  care  of 
through  modeling. 

A  state  space  realization  of  P  is  given  by 


A 

Sd 

B„ 

Ce 

Pad 

Do. 

Dcu 

Drd 

D.. 

Dru 

Py<i 

Pyw 

Pyu  . 

which  results  in: 


X  =  Aj,  +  *  B^w  *  B^u 

e  =  C^x  * 


(2.12) 


(2.13) 


(2.14) 
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(2.15) 


z  =  C,x  ^  ^  +  D,„u 

Y  =  C^x  +  Dyjd  +  +  D^^u 


Note  that  Dj,,  must  be  zero  or  the  2-norm  of  will  be 
infinite  for  any  controller.  and  are  assumed  to  be 
zero  for  convenience.  Further,  the  plant  must  satisfy  the 
following  conditions: 


1) 

(A,BJ 

stabilizable  and  (Cy,A)  detectable 

2) 

D 

D 

eu  ^eu 

full 

rank 

and  Dyd  full  rank 

3) 

D 

full 

rank 

and  Dy^,  DyJ  full  rank 

4) 

A-j,I 

Bu 

full  column  rank  for  all  w 

c. 

Beif 

5) 

A-j.I 

Bd 

full  row  rank  for  all  w 

c. 

Byd 

6) 

A-j«I 

Bu 

full  column  rank  for  all  w 

c. 

dJ 

7) 

j 

A-jwI 

B.| 

full  row  rank  for  all  w 

Cy 

DyJ 

Condition  1  must  be  satisfied  or  there  exists  no  stabilizing 
controllers.  Conditions  2)  through  7)  come  from  various 
aspects  of  the  individual  H,  and  H,  problems  which  are 
assumed  to  be  met  in  the  mixed  Hj/H.  problem  [1]. 
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The  compensator  K  is  given  in  state  space  by: 


x^=AcX^*B^y 


(2.17) 


u  =  C^x^ 


(2.18) 


Again,  to  avoid  an  infinite  must  equal  zero.  This 

does  not  result  in  loss  of  generality  [1].  This  thesis 
assumes  that  the  order  of  the  controller  equals  the  order  of 
the  plant.  It  as  been  shown  [4]  that  the  optimal  controller 
has  an  order  higher  than  that  of  the  plant  (possibly 
infinite).  It  is  often  true  that  higher  order  controllers 
gain  little  in  performance. 

Closing  the  P-K  loop  results  in: 

.  ,  .  -  (2.19) 

x  =  Ax  +  B^d*B^w 


(2.20) 


(2.21) 


where 
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Note  that  closed 
Now, 


where  is  the 


Also,  s  Y 


A  = 


B^Cy 


B.C, 

Ac 


B 


W 


B. 


C  =  \C  DC] 

<  '-e  ^ou^c  ‘ 

C  =  f  C  DC] 

‘-'zu^c  • 


(2.22) 

(2.23) 

(2.24) 

(2.25) 

(2.25) 


loop  system  stability  is  determined  by  A. 


Vl\X,  =  tr\Q,cJc^] 


(2.27) 


lolution  to  the  Lyapunov  equation 


(2.28) 


.f  there  exists  a  Q*  =  Q3  a  O  satisfying 

( 2 

AC?.  +  C?.A~^  +  Y  ^V^CjCaO.*B^Bj-:0 
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with  A  stable.  A  is  stable  if  (A,  B<i)  is  stabilizable . 

From  here,  the  problem  can  be  cast  as  a  Lagrange 
Multiplier  Problem  with  J  =  |Tj^l2  =  tr  [Q2^\C^'\  (minimizing 
the  square  of  the  function  is  the  same  as  minimizing  the 
function).  The  constraints  are  the  Lyapunov  Equation  2.28 
and  the  Ricatti  Equation  2.29. 

se=  tr  f  OaC/Cj  +  +  +  (2-30) 

^  tr{  [  AQ,  +  Q.A'^^y  -  B^J  ]  Y} 

The  first  order  necessary  conditions  for  a  minimum  are  found 
by  evaluating. 

=  =  =  =0  (2.31) 

dA^  dB^  dC^  5X  BY  dQ^  dQ. 

These  necessary  conditions  represent  a  set  of  7  matrix 
equations.  These  matrix  equations,  however,  do  not  lend 
themselves  to  an  immediate  solution.  Numerical  techniques 
must  be  used.  The  particular  technique  used  in  this  thesis 
is  covered  in  Chapter  3. 
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Chapter  III.  Analysis  and  Prcxredure 


3 . 1  Solution  Technique  [ 1 ] 

The  objective  is  to  generate  a  llT^^lj  versus  |T^|„  curve 
for  a  given  system.  This  entails  selecting  ay  (Yo  <  Y  < 

Y2)  and  solving  the  mixed  problem.  The  7  necessary 
conditions  mentioned  in  Chapter  2  are  coupled  and  highly 
nonlinear.  No  known  analytic  solution  exists,  so  a 
numerical  technique  must  be  used. 

In  order  to  facilitate  a  numerical  solution,  the 
performance  index  is  changed  to 

J=  (1-p)  trlO^cJc^]  IQ.CJC^]  (2.36) 

where  p  is  a  numerical  convergence  parameter.  At  p  =  1,  the 
problem  reverts  to  the  minimum  entropy  problem  with  a  well 
known  and  easily  obtainable  solution.  In  [1:134],  it  is 
shown  that  as  p-.  0  the  solution  does  converge  to  the  mixed 
Hj/H*  solution. 

The  choice  of  numerical  technique  is  critical,  since 
there  are  a  large  number  of  unknowns,  the  problem  is  highly 
nonlinear,  and  computation  of  second  derivatives  (fourth 
order  tensors)  is  unwieldly.  In  connectior.  with  research  in 
this  area,  a  Davidon-Fletcher-Powell  (DFP)  algorithm  had 
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already  been  coded  up  in  Fortran.  The  code  was  modified  to 
include  a  Ricatti  solver  and  to  allow  for  increased  problem 
size . 

The  DFP  problem  is  formulated  [  8  ]  as  minimizing  a  scalar 
function  F  given  a  vector  of  unknowns  X.  The  direction,  S, 
in  wnich  to  move  the  unknowns  is  given  as  a  function  of  the 
gradient  VF(x) 


S^H-vF{x) 


where  H  is  the  second  derivative  matrix.  H  does  not  have  to 
be  analytically  derived;  it  c~  .  -j  approximated  adequately. 
The  step  size  K  is  d  .termined  by  minimizing  F(X  +  KS)  by 
numerical  means.  Equality  boundaries  are  enforced  by 
artifically  setting  F(X  +  KS)  to  a  large  number  (10^*)  if 
the  boundary  is  violated.  This  has  the  effect  of  causing 
the  step  to  back  away  from  the  boundary. 

In  this  particular  problem,  the  boundary  constraints 
are  determined  by  four  of  the  seven  necessary  condition 
equations 
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AQ-  +  O^A  ■’  +  B^bJ  =  O 


(3.1) 


+  (1  -  Ji )  cjc^  =  0 

AQ^  +  O.aI  y  -'0-C/Ce 0-  +  5^/  =  O 
[  A  +  Y  ^Q.CjC^  ]  ^Y  +  Y  [  A  *y  ~^0.Cjc^  ]  +  ^  cjc^  =  O 


As  shown  in  [  1:14  3,144  ]  the  Q; ,  Q, ,  X  and  Y  solutions  must 
be  positive  semidef  inite .  Also,  for  fA#0 ,  (A  +  0.0,0^)  must 

be  stable.  Note  that  A  contains  the  controller  A^.,  and 
Cg  matrices,  and  thus  contains  all  the  unknowns.  The  number 
of  unknowns  is 

(#  states)^  +  (#  inputs )(#  states)  +  (#  outputs ) (#states ) 

Convergence  is  checked  by  observing  the  change  in  F  with 
respect  to  the  last  iteration's  F.  With  e  as  the 
convergence  criteria,  this  can  be  formulated  as 

vF  ^  (3.5) 

jF(x) 1 

The  flow  diagram  for  this  algorithm  is  given  in  Figure  3.1. 
This  problem  does  not  have  guaranteed  convexity,  so  a  good 
starting  point  is  critical.  The  starting  point  here  is  a 
guess  at  the  controller.  In  addition  to  satisfying  s 
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Y,  the  3tart  point  controller  must  satisfy  Equations  3.1 
through  3.4.  Typically,  when  starting  a  new  problem,  this 
would  be  chosen  to  be  the  H*  central  controller,  although 
any  admissible  "nearby"  controller  is  acceptable. 


19 


Figure  3.1  DFP  Flowchart 
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3.2  Procedure 


Matlab  is  used  to  initially  characterize  the  system  by 
finding  Yo»Y2»  *o*  Normally,  a  y  slightly  greater  than 

Yo  is  chosen,  and  then  Matlab  can  be  used  to  generate  the 
central  controller  as  a  start  point.  This  controller  is 
then  transferred  to  a  Fortran  input  file  for  use  by  DFP.  Of 
note  is  that  the  problem  is  easier  to  solve  at  y  close  to  Y2 
when  using  "overbound"  start  points,  because  the  start  point 
3c-norm  is  usually  far  less  than  the  y  that  generated  it. 

Near  Yo/  the  slack  in  the  |T^|sY  portion  of  the  problem  is 
getting  much  tighter  (in  fact  going  to  zero).  This  means 
that  the  start  point  becomes  increasingly  critical  as  Y-*Yo* 
At  each  Y  level,  the  p  that  is  initially  set  is  a  guess  that 
takes  into  account  how  close  the  starting  point  is  to  the 
minimum  entropy  solution  and  how  much  p  can  be  decreased 
from  the  starting  point  and  still  get  numerical  convergence. 
Starting  with  the  central  solution,  the  first  DFP  run 
would  probably  use  .5sps.99.  The  controller  output  from 
that  run  is  then  used  in  a  subsequent  DFP  run  at  a  lower  p. 
The  amount  that  p  can  be  decreased  between  runs  seems  to 
depend  largely  on  the  closeness  to  Yo*  Near  Yo»  might  be 
decremented  as  0.9,  0.7,  0.5,  0.4,  0.3,  0.2,  0.1,  0.01, 
whereas  near  Y2f  M  might  be  decremented  as  0.5,  0.1,  0.01, 
0.001.  This  is,  of  course,  problem  dependent,  and  one  must 
keep  in  mind  that  the  problem  may  be  changing  drastically 
with  changes  in  p.  From  the  problems  that  were  done  in 
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connection  with  this  thesis,  a  ^  of  about  0.001  is  adequate 
for  convergence  to  the  mixed  solution  at  a  given  y  level. 
This  gives  iT^dl.  =  y  to  about  four  decimal  places  and  iT^lj 
unchanging  to  about  7  decimal  places.  An  advantage  to 
starting  the  y  sweep  near  Yo  moving  toward  y,  is  that 
subsequent  y  levels  can  be  started  using  the  converged 
output  controller  from  the  previous  y  level  (iTedl^sy  and  K 
stabilizing  is  automatically  satisfied).  Thus,  the  H, 
central  controller  need  not  be  generated  for  each  y  level. 

A  procedural  flowchart  is  given  in  Figure  3.2. 
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Figure  3.2  Procedural  Flowchart 
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3 . 3  Set-u 


For  a  given  physical  system,  the  block  diagram  for  the 
system  is  manipulated  into  the  state  space  form; 


x^Ax  +  B^d  +  B^w  +  B^u 

e  =  C^x  +  D^d  + 

z  =  C^x  +  D^^d  + 

y=C^X  +  By^d  +  Dy^W  +  Dy^U 


There  are  several  tuning  parameters  available  for  the 
designer  to  control  the  optimization.  The  P  matrix  elements 
Bd/  Cg,  Cj,  Dyd/  Deu,  and  0,^  are  selectable  and  have 

the  following  interpretation: 
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A 

fid 

fi. 

fiu 

state  penalty-* 

Co 

0 

0 

Deu 

H2  state  penalty-* 

0 

0 

Dy^ 

0 

i  I 


03 

<0 

•■s 


<S) 

03 


O 

03 

C 


to 


*«i 

o 

03 


C 


03 


to 


control  penalty 
control  penalty 


From  here,  an  input  file  is  created  (see  Appendix  A  for  an 
excunple).  The  DFP  program  has  all  real  numbers  in  DOUBLE 
PRECISION  format.  For  this  thesis,  all  DFP  runs  were  made 
with  the  following  parameters: 


Q2  positive  definiteness  boundary  =  -10’^° 

Q,  positive  definiteness  boundary  = 

X  positive  definiteness  boundary  = 

Y  positive  definiteness  boundary  =  -10'*° 

(A  +  Y-2  Q®  Cg  Cg]  stability  boundary  =  +10‘® 
Convergence  Criteria  e  (Checkstop)  =  +10'® 
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Chapter  IV.  One  Exogenous  Input  One  Exogenous  Output 


4 . 1  Problem  Synthesis 

This  system  has  one  exogenous  input  and  one  exogenous 
output  (IIIO).  In  standard  form,  it  is  represented  as  in 
Figure  4.1: 


e 


y 


A 

p= 

Cc 

O 

Figure  4.1  II 10  Standard  Form 

Thus,  the  mixed  problem  is 

inf  |T^dl2  subj-ct  to  |T^d0-  ^  Y 

K  stabilizing 
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so  that  we  are  expecting  both  bounded  energy  inputs  and 
noise  inputs  to  be  contained  in  d  and  we  are  looking  at 
performance  and  robustness  with  respect  to  the  same  transfer 
function.  In  other  words,  performance  with  a  guaranteed 
level  of  robustness,  with  respect  to  the  same  transfer 
function . 

This  is  equivalent  to 

inf  ltT,„|lj  subject  to  «  Y 

K  stabilizing 

with  Ted  “  Tjy.  An  associated  block  diagram  for  such  a 
system  is  shown  in  Figure  4-2. 


Figure  4.2  System  Block  Diagram 
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4 . 2  Background 


Prior  to  the  Mixed  Solution,  the  best  solution  to  this 

problem  was  the  Minimum  Entropy  or  H*  central  solution.  In 

this  chapter,  the  Mixed  Solution  will  be  compared  with  an 

H*  central  solution  for  an  example  system.  It  has  been 

* 

shown  [6:12]  that  as  y  ->  *,  the  Minimum  Entropy  (ME)  and  the 
problem  are  equivalent.  For  y  less  than  «,  iTe^ll*  with 
the  ME  controller  is  less  than  y  and  |Tedl2  with  the  ME 
controller  is  an  upper  bound  on  the  minimum  This 

problem  has  an  analytic  solution.  Changing  the  block 
diagram  to  that  in  Figure  4.3,  K  is  the  lower  linear 
fractional  transformation  of  J  and  Q,  where 


J'= 

-^c 

O 

= 

5,  O  ^ 

J, 


uy 


oy 


J  is  found  by  solving  a  pair  of  algebraic  Ricatti 
equations  (with  both  solutions  positive  semidef inite ) .  The 
only  constraint  on  Q  is  that  it  be  stable  and  that  |Q|,  <  y. 
Thus,  a  K  such  that  iT^dl  ^  y  is  not  unique.  However,  in 
[6:19-22]  it  is  shown  that  the  ME  controller  is  just  J  with 
Q=0,  and  hence  the  name  H.  central.  The  ME  controller  then 
is 
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As  noted  in  Chapter  2,  this  K  is  the  same  as  the  Mixed 
Solution  with  n  =  1 .  Of  key  importance  is  that  the  ME 
solution  is  really  only  applicable  to  II 10  systems. 


Figure  4.3  Extended  Standard  Form 
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4 . 3  Description 

The  system  chosen  was  the  same  as  in  [ 1  ] .  It  is  a 
third  order,  unstable,  minimum  phase  system,  with 

Xi,  Xj  =  -.6812  ±  j  1.182 
X3  =  +.9549 

zeros  at  -1.1990 
-  .0490 

The  open  loop  singular  value  (SV)  plot  is  shown  in 
Figure  4.4.  Note  the  lightly  damped  zero.  The  state  space 
matrices  are: 


f-.3908  -.4565  1.266 
1.445  -1.049  1.208 
-.1288  .6744  1.032 


.  0488 

-.4275 

.3608 

.3564 

-.4470 

-.9172 

Cg  =  [.9420  .0144  .118/1 

Cy=[-1.557  -1.943  -.0914] 


.5185] 


1.357  5] 
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Magnitude 


10‘ 


Looking  at  the  closed  loop  system  using  the  optimal  H, 
controller,  it  can  be  seen  in  Figure  4.5  that  an  infinite 
bandwidth  system  results  with  =  2.1426.  As  shown  in 
Figure  4.6,  the  controller  also  has  infinite  bandwidth. 
Closing  the  loop  with  the  controller  results  in  a 

higher  oc-norm  and  a  more  typical  finite  bandwidth  singular 
value  plot  (Figure  4.5).  This  controller  results  in  y^  = 
4.163  and  =  3.248.  Figure  4.6  shows  that  this  controller 
has  finite  bandwidth.  Of  note  is  that  both  controllers  have 
a  spike  corresponding  to  the  lightly  damped  zero  in  the 
plant . 
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4 .4  Results 


The  Mixed  Solution  for  this  system  was  determined  using 
the  procedure  outlined  previously.  As  shown  in  Figure  4.7, 
the  Mixed  Solution  is  clearly  better  than  the  H,  central 
solution.  The  maximum  difference  between  the  two  solutions 
is  about  5  percent  at  a  of  about  2.5.  As  Figures  4.7 
and  4.8  show,  the  two  solutions  converge  as  y-^Yj  or 
Table  4.1  shows  the  comparison  for  selected  Y„i,  levels. 

Note  that  care  must  be  taken  in  such  a  table  because  y^i,  = 
llT^dL  but  |T^dll»  -  Yh  CENTRAL  SO  that  for  a  fair  comparison 
either  y.^,  or  Yh  central  should  be  adjusted  such  that  iT^dlo.  is 
the  same. 


Figure  4.7  Comparison  of  Mixed  and  Central  Controllers 
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I ' - 1 - ' - 1 - 1 - 1 _ I _ I 

2  2.2  2.4  2.6  2.8  3  3.2  3.4 

Ted  Infinity  Norm 

Figure  4.8  Mixed  and  Central  2-Nonn  Difference 


Table  4.1.  Comparison  of  Mixed  vs  Central 


Mixed 

H*  central 

Ymix 

1  2 

Y  BC 

1  Tgdl  • 

iTedh 

2.15 

2.15 

27.6012 

2 . 150 

2.15 

27.601 

2.5 

2.5 

4.3536 

2.540 

2.5 

4.5993 

2.8 

2.8 

3.6345 

2.9608 

2.8 

3.7956 

3.7 

3.7 

3.2594 

5.807 

3.7 

3.2765 

4.163 

4 . 163 

3.248 

00 

4.163 

3.248 

Complete  numerical  data  for  the  above  curves  can  be  found  in 
Appendix  B. 

Figure  4.9  shows  the  closed-loop  Mixed  Solution  results 
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for  Ymix  =  2.15,  2.5,  2.8,  and  3.7  as  well  as  the  and 

Hjopt  curves.  In  Figure  4.9  the  arrow  shows  decreasing  y. 
The  progression  from  H^opt  to  H^pt  with  decreasing  y  is 
clearly  evident.  Figure  4.10  shows  the  closed-loop  H, 
central  solution  results  following  the  same  trend  from  Hj^pt 
to  H,opf  Note  that  at  H^opt,  which  corresponds  to  y  =  the 
bandwidth  is  about  1  rad/s,  while  at  y  =  2.15  the  bandwidth 
is  about  300  rad/s.  The  bandwidth  goes  to  infinity  with 
just  a  small  change  to  y  =  2.1426,  illustrating  the  highly 
non-linear  nature  of  the  problem. 


Frequency  (rad/s) 


Figure  4.9  o-Plot  of  Mixed  Solution  T^^ 
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Frequency  (rad/s) 


Figure  4.10  o-Plot  of  H.  Central  Solution 

In  Figure  4.11/  the  Mixed  Solution  controllers  are 
shown  (solid  line  shows  and  for  the  same  y  levels 

as  before.  Again,  with  the  arrow  in  the  figure  showing 
decreasing  y  levels,  the  progression  from  Hj^pt  to  li^pt  is 
clearly  evident.  Figure  4.12  shows  very  simileur  results  for 
the  H.  central  controller. 


wagxiituuc  Magnitude 


Frequency  (rad/s) 

Figure  4.11  o-Plot  of  the  Mixed  Controllers 


Frequency  (rad/g) 

Figure  4.12  o-Plot  of  the  Central  Controllers 
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Comparing  the  Mixed  Solution  and  H,  central  solution  at 
specific  |T^|,  levels  reveals  some  of  the  differences 
between  the  two  solutions.  As  shown  in  Figures  4.13,  4.14, 
4.15,  and  4.16  there  is  virtually  no  difference  in  the 
controllers  below  1  rad/s.  It  is  at  the  higher  frequencies 
that  the  differences  become  apparent.  The  roll-off  appears 
to  be  similar  at  high  frequencies,  but  the  Mixed  Solution 
roll-off  appears  offset  downward.  The  amount  of  off-set  is 
greatest  in  the  iTe^l.  =  2.5  to  2.8  range,  which  corresponds 
to  the  area  of  greatest  IT^I^  difference  between  the  two 
solutions.  The  solutions  become  very  similar  at  y  near  y^ 
and  y  near  y^,  as  expected. 


Fiequency  (rad/s) 


Figure  4.13  o-Plot  of  Controllers  (oc-norm  =  3.7) 


38 


Frequency  (rad/s) 


Looking  at  the  closed  loop  comparisons  in  Figures  4.17, 
4.18,  4.19,  and  4.20,  we  see  a  very  interesting  phenomena. 

At  low  frequency,  the  curves  are  flat  (at  the  value), 

and  then  a  roll-off  begins  near  1  rad/s.  At  this  roll-off, 
however,  the  Mixed  Solution  appears  to  make  a  sharper  corner 
than  the  H,  central  solution  and  its  singular  value  plot 
then  passes  under  the  H*  central  curve.  This  characteristic 
is  especially  evident  in  Figures  4.18  and  4.19.  Thus,  the 
mixed  solution  has  a  lower  2-norm  because  it  is  trading  low 
frequency  gain  for  high  frequency  gain  to  lower  the  overall 
area  under  the  curve.  The  Mixed  Solution  is  able  to  do  this 
because  it  can  make  a  sharper  corner  to  the  roll-off.  In 


40 


Magnitude 


fact,  in  [4],  it  is  conjectured  that  the  optimal  order 
controller  would  have  a  discontinuous  corner.  This  would 
take  an  infinite  order  controller. 


Frequency  (rad/s) 

Figure  4.17  o-Plot  of  ^x-norm  =  3.7) 


4 


Magnitude  Magnitude 


Frequency  (rad/s) 

Figure  4.18  o-Plot  of  T«,  (oo-norm  =  2.8) 


Frequency  (rad/s) 

Figure  4.19  o-Plot  of  (oo-norm  =  2.5) 


Magnitude 


Frequency  (rad/s) 

Figure  4.20  o-Plot  of  (oo-norm  =  2.15) 


Thus  for  the  II 10  system  the  Mixed  Solution  is  superior 
to  the  central  solution.  As  expected,  the  two  solutions 
converge  at  and  Y2*  The  difference  between  the  two 
solutions  seems  to  manifest  itself  as  a  sharper  corner  to 
the  roll-off  in  the  o-plot.  This  difference  is  greatest 
for  y's  roughly  mid-way  between  Yo  Yz*  The  next  level 
of  problem  is  the  1120  system  that  will  be  covered  in 
Chapter  5 . 
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Chapter  V.  One  Exogenous  Input  Two  Exogenous  Outouts 


5 . 1  Problem  Synthesis 

This  system  has  one  exogenous  input  and  twc-  axogenous 
outputs  (1120).  In  standard  form  it  is  represented  as  seen 
in  Figure  5.1: 


e 

z 


y 


Figure  5.1  1120  Standard  Form 

Thus ,  the  mixed  problem  is 

inf  I|T,„D2  subject  to  IT^I.  Y 

K  stabilizing 

so  that  we  are  modelling  both  bounded  energy  inputs  and 
noise  inputs  in  w  while  evaluating  the  performance  and 
robustness  with  respect  to  separate  outputs.  A  practical 
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motivation  for  such  a  problem  might  be  an  active  control 
surface  on  an  airplane;  the  motion  of  the  airplane  would 
respond  to  actuator  noise  (or  perhaps  stabilizing  high 
frequency  commands)  as  well  as  typical  step  flight  path 
commands  coming  into  the  control  surface.  Then  and  Te„ 
might  represent  high  and  low  frequency  models,  respectively. 

5 . 2  Background 

The  Mixed  Solution  has  full  capability  to  handle  a  wide 
range  of  different  models  resulting  in  a  1120  system. 
However,  prior  to  the  Mixed  Solution  the  only  solution  to 
this  problem  was  Bernstein  and  Haddad's  (B  &  H)  approach 
forwarded  in  [5].  This  approach  uses  an  overbound  technique 
and  is  somewhat  restrictive,  as  will  become  apparent.  The 
Mixed  Solution  will  be  compared  with  a  B  &  H  solution  on  an 
example  system  shown  in  [5].  The  general  B  &  H  block 
diagram  is  shown  in  Figure  5.2. 
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62  Zj  Wj  Cl  Z, 


Figure  5.2  System  Block  Diagram 

Briefly  reviewing  the  B  &  H  approach,  the  equivalent 
B  &  H  state  space  representation  is 


A 

E, 

B 

E,. 

O 

E2. 

E, 

o 

E2 

C 

0 

and  the  following  coupled  set  of  Ricatti  equation  must  be 
solved ; 
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O^AQ^OA  ^  +  Q* 


(5.1) 


0=  (A  +  y  Mo  +  i?]  i?i-)  ^P+P(A  +  y  MC)  +  £?]  (5.2) 

-5^PEPS  +  J?j 

0=  (A-ZP  +  y  ^OP,-)  (?+(?(A-ZP  +  y  (5.3) 

+  (?(Y  ^  (Pi.  +  P^S^PSPS)  )  C+C^Q 


where 


y,  -AD/ 

r,  =  e/e^ 

£  =  c^(D2A^)  c 

a.=pa 

5=  (i  +  p^Y  2£^P)  1 


Note  that  the  H*  control  penalty  must  be  a  scalar 
multiple  of  the  Hj  control  penalty.  The  general  B  &  H 
equations  5.1,  5.2,  and  5.3  are  not  analytically  solvable, 
and  thus  far  no  one  has  numerically  solved  them  either. 
However,  with  p  =  0  they  are  solvable  by  Homotopy 
Techniques .  This  is  the  case  where  there  is  no  control 
penalty,  and  in  the  example  problem  presented  in  [5]  this  is 
what  was  chosen.  The  Ricatti  Equations  then  become: 
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0  =  AC>  +  C>A’‘  +  C>(Y“^/2i.-£)  0+ V^i 

o-  (A*y^[Q*C]  R^lP  +  PiA^y-^  10^6)  R^J 

~PEP+R^ 


(5.4) 


(5.5) 


0=  (A-SP  +  y  6*(5{A-LP*y^QR^J'^  (5.6) 

+  i?(Y  (?+(?£  P 


where 


c^=-  (P2) 

A^=A  +  BC^-flpC  +  Y  ^P^i- 


Eqn  5.4  is  uncoupled  and  may  be  solved  immediately.  A 
direct  Homotopy  approach  will  solve  Equations  5.5  and  5.6 
for  y*  <  y  <  <x  where  Yo  <  Y*  Y2r  y*  may  be 

substantially  larger  than  y^.  Therefore,  a  full  vs 

iTewIx  plot  cannot  be  generated.  A  direct  Homotopy  algorithm 
was  coded  in  Matlab  M.file  form  (BH.M  in  Appendix  A),  but  it 
was  not  possible  to  get  y*  very  low  for  any  example  tried. 

A  different  algorithm,  forwarded  in  [3],  was 
subsequently  coded  up  ( BLACK. M  in  Appendix  A)  and  found  to 
be  much  better.  The  Ricatti  equations  for  this  algorithm 
are: 
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0=A-^P^-^P^A-P^LP^*R^ 

0=  (A-SP)  (ASP)  ^*0^0 

0  =  AQ^QA'^*Q(r^R^^-2)  0+ Vj 
0=  (A  +  y  "  tO  +  C?]  ^,J  ’■P  +  PU  +  Y'^  [O^C?]  PjJ 

-PSP  +  Pj 

0=  (A-EP  +  y  -y)  ^1-)  1 

(a-Ep  +  y‘M£?+ 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 


The  flowchart  diagram  for  this  algorithm  is  shown  in 
Figure  5.3.  Of  note  is  that  the  only  input  required  in  this 
algorithm  is  y.  The  algorithm  generates  its  own  start  point, 
so  "solution  cascading"  going  from  a  higher  to  a  lower  y  is 
not  required.  This  algorithm  was  used  to  generate  all  B  &  H 
solutions  in  this  chapter. 
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findP„(5.7) 


Figure  5.3  Black  Algorithm 
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5.3  The  Dual 


For  completeness,  the  dual  of  this  problem  must  be 
covered.  This  problem,  as  represented  in  Figure  5.4,  was 
first  forwarded  in  [7].  It  was  subsequently  shown  in  [2]  to 
be  the  dual  of  the  BH  problem. 


Figure  5.4  Dual  Standard  Form 

This  problem  is  formulated  as 

inf  iTpwlj  subject  to  |T^|,  s  y 

K  stabilizing 

so  that  we  are  modelling  bounded  energy  inputs  in  d  and 
noise  inputs  in  w,  with  all  outputs  modelled  in  e. 

The  mixed  state  space  representation  is 


A 

Ba 

B^ 

B. 

p= 

Co 

O 

0 

Beu 

c. 

Dy. 

By. 

0 
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while  the  equivalent  B  4  H  state  space  is 


A 

Pi.  Pi 

B 

p= 

Er 

O  0 

E2 

C 

P2.  P2 

0 

Similar  to  the  original  B  &  H  problem,  must  be  set  to 
zero  to  obtain  the  following  set  of  Ricatti  equations; 


0=  {A  +  y-^V^^lO*d])  P*P{A*y-^V^^[0*d]  )  ^  (5.13) 

+  p2p+ 


0=  (A-pS+y’^ViJ:?)^Q  +  P(A-p2+y-^Vi,C>)  (5.14) 


where  all  definitions  are  the  same  as  before,  but  with  the 
addition  of  Dj*,  and  Vj,  (Vj*  =  Dj:.  DjI).  A  direct 

Homotopy  algorithm  in  Matlab  M.file  form  (BHD.N  in  Appendix 
A)  was  coded,  but  no  comparisons  with  the  Mixed  Solution 
were  done  for  this  dual  problem. 

5 . 4  Description 

The  system  chosen  was  the  same  as  in  [5].  It  is  an 
eighth  order,  neutrally  stable,  non-minimum  phase  system, 
with ; 
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Xi,  =  -.037  ±  j  1.85 

Xj,  =  -.028  ±  j  1.41 

K,  =  -.015  ±  j  7.65 

K,  K  ^  0,  0 

zero  at  1.0 

In  [5],  it  is  said  that  this  represents  a  physical  system 
which  has  coupled  rotating  disks  with  noncolocated  sensors 
and  actuators. 

The  open  loop  singular  value  plot  is  shown  in  Figure 
5.5.  Of  special  note  is  the  response  between  0.5  and  3 
rad/s,  with  three  ripple-like  spikes  at  0.7,  1.5,  and  1.9 
rad/s.  As  will  be  seen,  the  controller  is  most  active 
controlling  these  spikes.  The  state  space  matrices  are: 

-.161  1000000 
-6.004  0100000 
-.5822  0010000 

-9.9835  0001000 
-.40/3  0000100 

-3.982  0000010 

0  0000001 

0  0000000 


0000  .00055  .011  .00132  .018 

0000  0  0  0  0 

0000  .00055  .011  .00132  .018 

0000  0  0  0  0 

Cy  =  C=  [100000001 

=  fO  1  ] 
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0 

0 


0 


0 

.0064  0 
.00235  0 
.0713  0 
1.0002  0 
.1045  0 


0 

0 

.0064 

.00235 

.0713 

1.0002 

.1045 

.9955 


.9955  0 


Looking  at  the  form  of  the  above  matrices  we  see  that 
an  equivalent  block  diagram  would  be  like  that  in  Figure 
5.6.  The  Wj  disturbance  is  coming  into  the  plant  added  to 
u.  The  disturbance  Wj  is  coming  into  the  output  y. 

Finally,  the  outputs  ej  and  Zj  are  being  collected  (through 
the  weighting  matrix  p„)  directly  from  the  plant  states. 
Thus  the  plant  output  y  is  uncoupled  from  the  "unmeasured" 
ej  and  Zj  outputs  that  have  oo-norm  and  2-norm  significance, 
respectively.  Note  the  absence  of  e,  to  place  a  penalty  on 
control  effort  expended  on  the  oo-norm  side  of  the  problem. 
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ISO 


Because  the  H,  control  penalty  is  zero,  the  HL  optimal 
controller  is  difficult  to  calculate  for  this  problem.  In 
[3],  a  Yo  of  0.12  was  found  for  this  problem  using  bisection 
techniques.  Closing  the  loop  with  the  controller  gives 
a  with  the  characteristics  shown  in  Figure  5.7.  This 
controller  results  in  =  0.37856.  Note  that  this 
controller  (Figure  5.8)  has  three  distinct  dips 
corresponding  to  the  0.7,  1.5,  and  1.9  rad/s  ripples  in  the 
upper  o-plot.  Figure  5.9  shows  with  this 
controller.  Here,  Y2  =  |Tewl-  =  1.3923.  Because  ej  and 
share  the  same  weighting  matrix  (p„  in  Figure  5.5),  the 
lower  o-plot  is  very  similar  to  the  o-plot  with  the 
same  spikes  and  dips. 
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Magnitude  Magnitude 


Frequency  (rad/s) 

Figure  5.8  o-Plot  of  the  Hjopt  Controller 


Frequency  (rad/s) 


Figure  5.9  o-Plot  of  the  with  T„  Controller 
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Since  the  H„  central  controller  is  diffir-ult  to 
generate  for  this  problem,  the  B  &  H  results  served  as  a 
starting  point  for  the  Mixed  Solution.  This  worked  very 
well  as  long  as  Ybix  significantly  greater  than  |Te„|„BH. 
The  Mixed  Solution  would  reliably  start  with  fx  =  0.01  in 
this  case.  The  problem  is  that  when  the  B  &  H  problem 
cannot  be  solved  at  a  y  near  the  Mixed  Solution  has  no 
start  point.  The  practical  limit  of  about  y  =  0.27  was 
found,  below  which  the  B  &  H  solution  wouldn't  work  to  start 
the  Mixed  Solution. 

5.5  Results 

The  Mixed  Solution  is  clearly  better  than  the  B  &  H 
solution  as  shown  in  Figure  5.10.  As  y^  is  approached,  the 
difference  between  the  two  solutions  becomes  large,  while  as 
y^  is  approached,  the  two  converge  to  the  same  answer. 

Table  5.1  shows  the  comparison  for  selected  y^j*  levels. 

Once  again,  care  must  be  taken  in  such  a  table  because  |Te»,|„ 
s  Ybh  while  |Te„U  =  Ymixf  so  either  or  y„i^  must  be  adjusted 
such  that  |Te„|*  is  the  same. 
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Tew  Infinity  Norm 

Figure  5.10  Comparison  of  Mixed  and  B  &  H  Controllers 


Table  5.1  Comparison  of  Mixed  vs  B  &  H 


Mixed 

B  &  H 

Y  mix 

iTeJ. 

ilT^h 

Y  BB 

mam 

1  T,  J  2 

.27 

.27 

.8382 

.2750 

.27 

1.1531 

.40 

.40 

.5444 

.4133 

.40 

.61521 

.65 

.65 

.4281 

.7046 

.65 

.4554 

.75 

.75 

.4097 

.8412 

.75 

.4288 

1.00 

1.00 

.3865 

1.329 

1.00 

.3931 

1.37 

1.37 

.3786 

6.81 

1.37 

.3786 

1.3923 

1.3923 

.3786 

00 

1.3923 

.3786 

Complete  numerical  data  for  the  B  &  H  curve  can  be  found  in 
Appendix  B. 

Figure  5.11  shows  the  closed-loop  Mixed  Solution 
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results  for  Ymi*  =  0.27,  0.40,  0.65,  0.75,  1.00,  and  1.37  as 
well  as  the  H^opt  curve.  In  the  figure,  the  arrow  shows 
decreasing  y  •  The  progression  away  from  is  readily 

apparent.  Note  that  as  y  progresses  to  lower  values,  the 
major  change  in  is  at  low  frequency  (below  0.8  rad/s), 
and  above  1.5  rad/s  no  change  is  discernable.  In  Figure 
5.12,  the  closed-loop  Mixed  Solution  results  for  the 
same  y  levels  are  shown.  Again  the  smooth  progression  away 
from  Hjopt  is  evident.  Note  how  low  frequency  energy  is 
being  pushed  out  to  higher  frequencies  as  y  is  decreased. 

In  Figure  5.13,  the  Mixed  Solution  controllers  are  shown  for 
the  same  y  levels.  Of  note  here  is  that  the  Mixed 
controller  has  only  two  dips  in  the  y  =  1.00,  0.75,  and  0.65 
cases;  the  1.9  rad/s  dip  disappears. 
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Magnitude 


Frequency  (rad/s) 

Figure  5.11  o-Plot  of  the  Mixed  Te„ 


Frequency  (rad/s) 


Figure  5.12  o-Plot  of  the  Mixed 


Figures  5.14  and  5.15  show  the  closed-loop  B  &  H 
solution  and  results  for  the  same  y  levels.  These 
curves  show  the  seime  progression  away  from  Hjopt  with 
decreasing  y,  and  show  similar  trend  characteristics  as 
compared  with  the  Mixed  results.  Figure  5.16  shows  the  B  & 
H  controllers  for  these  same  y  levels.  It  is  very 
interesting  that  the  B  4  H  controllers  always  possess  three 
dips,  very  similar  to  the  controller  while  the  Mixed 

controller,  transitions  to  and  from  a  three  dip  controller 
depending  on  the  y . 
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Frequency  (rad/s) 

Figure  5.16  o-Plot  of  B  &  H  Controllers 


Comparing  the  Mixed  Solution  and  B  &  H  solution  at 
specific  llTg^lL  levels  reveals  some  of  the  differences 
between  the  two  solutions.  Starting  with  Figures  5.17  and 
5.18,  which  show  Te„  and  T„,  respectively,  at  [|T„|(,  =  1.37, 
we  see  the  two  solutions  are  very  close.  Although  it  can 
not  be  seen  on  these  two  plots,  the  Mixed  has  a  slightly 
greater  bandwidth  than  the  B  &  H  Te„,  while  the  Mixed  has 
a  slightly  higher  peak  sigma  value  than  the  B  &  H  T^^,.  Both 
of  these  trends  will  continue  as  y  is  decreased.  The 
controllers  shown  in  Figure  5.19  are  almost  identical  and 
both  show  the  three  dip  phenomenon.  Moving  down  to  II  T^^IL  = 
1.00,  Figure  5.2  0  shows  how  the  Mixed  has  a  larger 
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bandwidth  than  the  B  fc  H  T„.  As  in  Chapter  4,  it  will 
become  increasingly  apparent  that  as  y  decreases,  as  the 
Mixed  Solution  is  making  a  "sharper  corner"  to  roll-off  for 
the  00-norm  side  of  the  problem.  In  this  figure,  also  note 
the  spike  at  0.7  rad/s.  As  y  gets  lower  this  spike  will 
become  increasingly  important.  Figure  5.21  shows  that  the 
Mixed  Solution  is  reducing  higher  frequency  response  of 
by  increasing  lower  frequency  response  (at  the  0.08  rad/s 
hump).  Note  that  the  Mixed  controller  in  Figure  5.22  has 
transitioned  to  a  two  dip  controller,  while  the  B  &  H 
controller  still  has  an  Hjop^-like  three  dip  controller. 

Thus,  the  Mixed  Solution  is  not  dealing  with  the  1.9  rad/s 
ripple.  In  fact,  the  B  &  H  controller  always  appears  to 
!e  a  scaled  controller.  At  |Te„|*  =  0.65,  Figure  5.23 

shows  the  increasingly  large  spike  at  0.7  rad/s.  Note 
that  above  0.7  rad/s  both  the  B  &  H  solution  and  the  Mixed 
I'olution  are  the  same.  Figure  5.24  shows  the  increasing 
livergence  between  the  two  solutions  for  between  0.03 
.nd  10  rad/s.  Figure  5.25  shows  the  controller  comparison 
•^iimilar  to  the  |Te„|*  =  1.00  controller.  For  =  0.4 

figures  5.26  and  5.27  show  Te„  and  respectively.  Note 
that  for  the  0.7  rad/s  spike  nearly  reaches  |Te„|«.  Again 
note  the  sharpness  of  the  roll-off  corner.  In  the  curve 
we  see  the  Mixed  Solutions  pronounced  ability  to  reduce  the 
gain  above  0.3  rad/s.  This  was  achieved  with  only  a  small 
increase  in  the  maximum  T^^  gain.  Figure  5.28  shows  the 
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controller  comparison  for  the  |T„S«  level.  Interestingly, 
the  Mixed  controller  has  now  transitioned  back  to  three 
dips,  probably  to  deal  with  the  emerging  1.9  rad/s 
spike.  At  |Te„IL  =  0.27,  the  lowest  level,  the  problem  is 
changing  radically.  In  Figure  5.29,  the  0.7  rad/s  spike 
has  joined  with  the  Mixed  curve  at  the  roll-off  point 
resulting  in  a  very  sharp  roll-off  corner.  Meanwhile,  the  B 
&  H  roll-off  corner  is  quite  rounded.  For  the  first  time, 
higher  frequency  differences  are  apparent  as  well.  This  is 
shown  better  in  Figure  5.30.  The  T„  curve  in  Figure  5.31 
tells  us  why  the  I  T^^ll  j  are  so  different  for  these  solutions 
at  this  |Te,,U  level;  the  B  &  H  T„  has  a  large  gain  step  at 
0.7  rad/s  that  the  Mixed  Solution  just  does  not  have.  The  B 

&  H  solution  seems  unable  to  trade-off  this  large  spike  for 

a  slight  increase  in  lower  frequency  gain.  Figure  5.32 
shows  both  controllers  having  the  same  three  dip 
characteristic  as  Hjop,..  Thus,  we  can  conclude  that  the  Hj 
overbound  used  in  the  B  &  H  solution  becomes  increasingly 
poor  as  Y  is  decreased.  This  is  especially  evidenced  by  the 

llTewll.=  .27  T„  plot,  where  the  Mixed  Solution  is  making  shrewd 

high/low  gain  trade-offs  and  the  B  &  H  solution  is  not. 


66 


Magnitude 


(solid=B&H  dashed=Mixed) 


10-51— 

10-3 


I  1  lij-u _ t \ — I  1  HIM  -L.,,  I — »  ■  I  »  •  I  nn  _L_1_] 

10-2  10-‘  10°  IQi 

Frequency  (rad/s) 

Figure  5.23  o-Plot  of  (oo-norm 


Figure  5.24  o-Plot  of  (oc-norm 


Magnitude  Magnitude 


Magnitude  Magnitude 


Frequency  (rad/s) 

Figure  5.29  o-Plot  of  Te„  (oo-norm 


10'3  10-2  10-‘  100  IQl 


Frequency  (raH/s) 

Figure  5.30  o-Plot  of  (ac-norm 


Figure  5.32  o-Plot  of  Controllers  (sc-norm  =  0.27) 
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Chapter  VI.  Two  Exogenous  Inputs  Two  Exogenous  Outputs 


6 . 1  Problem  Synthesis 

This  system  has  two  exogenous  inputs  and  two  exogenous 
outputs  (2120).  In  standard  form,  it  is  represented  as  in 
Figure  6.1. 


e 

z 


P  = 


A 

Bu 

Ce 

0 

Be. 

Beu 

0 

B,. 

Byv 

0 

Figure  6.1  2120  Standard  Form 

Thu.:,  the  mixed  problem  is 

inf  iT^wh  subject  to  iT^dll*  ^  Y 

K  stabilizing 

so  that  we  are  optimizing  noise  performance  in  the  T„ 
transfer  function  while  we  are  concerned  about  robustness  in 
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for  example,  could  be 


the  transfer  function.  T^, 
associated  with  connnand  following  error,  while  could  be 
associated  with  plant  response  to  noise. 

6 .2  Background 

There  are  currently  no  other  techniques  to  handle  a 
2120  system.  The  Mixed  Solution  is  the  only  real  solution 
in  this  case.  It  is  possible  to  generate  an  H,  central 
solution  based  only  on  Tg^r  but  this  is  not  really  mixed  at 
all  since  Tj„  is  not  considered  in  the  problem.  The  only 
true  comparisons  in  this  case  are  with  the  and 

solutions  at  each  end  of  the  trade-off. 

6 . 3  Description 

The  B  &  H  system  in  Chapter  5  can  be  readily  extended 
to  a  2120  problem  with  a  block  diagram  as  in  Figure  6.2. 
Manipulating  the  state  space  we  have 

x  =  Ax  +  B(u  +  prfj  d, + 


^2. 

Pe2^. 

PzJ-* 

P^2^ 

or 
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flpd,  0  Bp^,  0  B 


PeJ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Pe2 

9mi 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

9^2 

C 

0 

Pd2 

0 

Pw2 

0 

The  A,  B„,  By,  C^, 
Chapter  5,  while 


Cjf  Cy,  Dy„,  and  D^y  matrices  are  as  in 


^eu- 


0 

.  0001 


^yd  “  ^yw 


Thus,  we  have  disturbances  entering  the  system  at  the 
plant  input  and  output.  Robustness  is  being  measured  with 
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respect  to  some  plant  states  different  than  the  measured 
states.  Noise  is  also  coming  into  the  system  at  the  plant 
input  and  output,  with  performance  being  measured  with 
respect  to  some  plant  states  different  than  the  measured 
states.  All  oi'  the  p  weighting  matrices  are  non-zero  in 
this  example. 

The  and  Hj^pt  controllers  are  as  shown  in  Figure 

6.3.  Note  the  large  gain  increase  in  the  controller 

above  1  rad/s  in  contrast  to  the  Hjop^'s  roll-off.  These 
controllers  result  in  Yo  =  0.12015,  y?  =  1.3923,  and  oc^  = 
0.3786.  Figure  6.4  shows  the  resulting  closed  loop  T^„  and 
Ted  with  the  Hjopt  controller.  Since  this  is  the  T^^-H^opt 
rontroller,  iT^^lj  is  minimized  here,  but  Ted  is  showing 
similar  characteristics  with  this  controller.  Figure  6.5 
shows  the  resulting  closed  loop  Ted  with  the  H^p^ 

controller.  Here  we  see  the  typical  flat  o  (Ted)  curve. 
Note  that  |Tj„|j  is  going  to  infinity. 

6 . 4  Results 

The  iTedI,  versus  IT^I^  trade-off  is  clearly  shown  in 
Figure  6.6.  This  data  is  also  displayed  in  Table  6.1. 
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Magnitude  Magnitude 


Table  6.1  Mixed  Solution  to  2120 


Y 

iTedl. 

|T„l2 

.27 

.27 

.8381 

.4 

.4 

.5444 

.7 

.7 

.4180 

1.1 

1.1 

.3824 

Figure  6.7  shows  the  controllers  transitioning  from 
Hjopt  to  with  the  arrow  pointing  toward  decreasing  y  • 

Note  the  large  change  from  y  =  .27  to  y^.  In  Figure  6.8,  Ted 
is  shown  progressing  to  the  infinite  bandwidth  solution 

as  y  is  decreased.  Note  that  while  |Tedl®  is  going  down,  the 
minimum  singular  value  is  going  up.  Thus,  a  trade-off  is 
being  made  here,  as  well  as  between  the  2-nonn  and  oc-norm. 
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In  Figure  6.9,  is  shown  progressing  away  from  the  H 
solution  as  y  is  decreased.  Note  that  the  bandwidth 
increases  here  as  y  is  decreased. 


Magnitude 


Looking  at  the  y  =  *27  case.  Figure  6.10  shows  and 
the  Mixed  controller  (shifted  up  two  decades  for  clarity). 
The  controller  peaks  at  .9  and  1.8  rad/s  can  be  seen  to 
correspond  with  dips  in  the  maximum  singular  value  at 
those  frequencies.  Notice  the  similarity  the  controller  has 
to  the  minimum  Tej  singular  value  plot,  caused  in  this 
problem  by  the  way  the  state  and  control  penalties  are 
structured.  Figure  6.11  shows  and  the  Mixed  controller, 
but  the  relationships  here  are  not  so  clear.  Figure  6.12 
and  6.13  show  similar  characteristics  for  y  =  .7. 


10-3  10-2  10-1  lOO  10*  102  103 

Frequency  (rad/s) 

Figure  6.10  o-Plot  of  Controller  and  (oo-norm  =  .27) 
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S  10-2 


Frequency  (rad/s) 

Figure  6.11  o-Plot  of  Controller  and  T,„  (»_norin 


Thus,  the  Mixed  Solution  provides  complete  visibility 
into  the  oo-norm  versus  2-norm  trade-off  for  the  2120  system 
This  practical  example  again  showed  the  solution  to  be  on 
the  oo-norm  boundary  where  8 Ted!*  =  Y* 
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Chapter  VII.  Conclusions  and  Recommendations 


7 . 1  Summary  and  Conclusions 

The  objective  of  this  thesis  was  to  apply  a  newly 
developed  nonconservative  mixed  optimization  solution  to 
several  excunple  problems  and  compare  this  solution  to  the 
best  previously  available  solutions.  The  final  measure  in 
the  comparison  was  the  2 -norm  versus  oc-norm  trade-off  curve. 

In  Chapter  2,  various  aspects  of  control  theory  and  the 
role  of  H,  and  H*  theory  were  discussed.  Chapter  3  covered 
the  procedure  and  analysis  techniques  required  to  obtain  a 
Mixed  Solution. 

The  1110  example  was  looked  at  in  Chapter  4.  The 
comparison  in  this  case  was  done  between  the  Mixed  and 
H.,  central  solutions  for  a  third  order  SISO  system.  The 
example  showed  that  the  Mixed  Solution  does  converge  to  the 
Hjopt  solution  at  Yj  and  to  the  solution  at  y^.  For  Yo  < 

Y  <  Y2f  the  Mixed  Solution  is  superior  to  the  H,  central 
Solution  because  of  the  latter's  conservative  overbound.  In 
this  specific  example,  the  difference  peaked  at  the  mid¬ 
range  Yf  and  was  approximately  5  percent  of  the  2-norm. 

A  1120  example  was  covered  in  Chapter  5.  For  this 
system,  the  Mixed  Solution  was  compared  with  the  Bernstein  & 
Haddad  Solution  for  an  eighth  order  system.  This  particular 
example  problem  came  from  the  original  Bernstein  &  Haddad 
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paper  demonstrating  their  solution  and  included  the 
limitation  that  the  H*  control  penalty  must  be  zero.  The 
example  showed  that  both  the  Mixed  Solution  and  the 
Bernstein  &  Haddad  solution  converge  to  the  same  (Hjop^) 
solution  at  Yj.  For  Y  ^  Y2/  the  Mixed  Solution  is  superior 
to  the  Bernstein  &  Haddad  solution,  again  because  of  the 
latter's  conservative  overbound.  As  in  the  1110  comparison, 
the  hallmark  of  the  Mixed  Solution  appears  to  be  the  sharper 
corner  to  the  roll-off  that  it  makes  on  the  plot. 

Chapter  6  looked  at  a  2120  eighth  order  example  system. 
This  system  was  chosen  as  an  extension  of  the  1120  Bernstein 
&  Haddad  system  because  of  its  practical  motivation 
(actuators,  sensors,  and  spinning  masses).  No  comparison  to 
previous  results  was  possible  in  this  case,  because  no  other 
solution  exists  to  this  problem.  The  example  again  showed 
again  that  the  Mixed  Solution  converges  to  the  H,  ^p.  or  H^ 
solution  at  their  respective  y's.  The  example  clearly 
showed  the  00-norm  versus  2-norm  trade  that  is  possible. 

7 .2  Recommendations 

Practical  applications  still  need  to  be  looked  at  in 
earnest.  Several  of  the  practical  examples  that  were  looked 
at  could  not  be  done  with  the  software  code  used  in  this 
thesis,  because  either  the  problem  would  not  numerically 
start,  or  once  started  would  not  converge.  It  is  not  known 
whether  this  is  related  to  modeling,  numerics,  or  software, 
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but  this  is  definitely  an  area  of  further  study. 

As  an  adjunct  to  this,  further  refinements  to  the 
current  numerical  algorithm  might  help  its  robustness  and 
speed.  Perhaps  a  whole  new  algorithm  could  better  find  the 
solution. 

For  the  current  algorithm,  finding  a  start  point  for  a 
system  greater  than  II 10  is  very  time  consuming  because  of 
numerous  points  that  should  work,  but  do  not.  A  reliable 
method  of  finding  start  points  would  be  a  worthwhile 
investment . 

Finally,  the  Mixed  Solution  has  a  large  number  of 
tunable  parameters.  The  practical  techniques  of  using  these 
parameters  to  achieve  desired  results  still  needs  to  be 
investigated . 
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Appendix  A:  Selected  Software  Items 

This  appendix  contains  the  PRO-MATLAB  user  written 
routines  utilized  in  this  thesis  and  an  example  of  an 
input  file  to  the  Fortran  based  DFP  program. 


BH.M  algor ithvim  for  solving  the  Bernstein  and  Haddad 
problem: 

g2=l/( gara*gam) ; 

a=am; 

b=bm; 

c=cm; 

vl=di*dl ' ; 
v2=d2*d2 ' ; 
rl=el ' *el ; 
r2  =  e2 ' *e2  ; 
rlinf=elinf'*elinf; 
sig=b*inv{ r2 ) *b' ; 
sigbar=c' *inv( v2)*c; 

% 

%  solve  for  q 
% 

f=a'  ; 

g=-( g2*rlinf-sigbar ) ; 
h=vi ; 

q=are( f ,g,h) ; 

% 

%initialize  p 
% 

%p=eye(n) ; 
p=pchest ; 

% 

%  iterate  between  qhat  and  p  equations 
% 

er  ror==l ; 

while  error>.001 
% 

f=(a+g2*q*rlinf-sig*p) ' ; 
g=-g2*rlinf ; 
h=q*sigbar*q ; 
qhat=are( f ,g,h) ; 

% 

f=(a+g2*{q+qhat)*rlinf ) ; 

g=sig; 

h=rl ; 

pold=p ; 

p=are( f ,g,h) ; 
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% 

error»norm( p-pold, 'fro'  ) 
% 

end 

disp( 'converged' ) 

% 

%  check  answers 
% 

if  abs(q'-q)>10“(-5) 
disp('q  not  symetric') 
end 


if  rain  (eig(q))  <  0 
disp('q  not  positive  definite') 
min(eig(q) ) 
end 
% 

if  abs(qhat'-qhat)>10"(-5) 
disp('qhat  not  positve  definite') 
end 
% 

if  min( eig ( qhat ) ) <0 
disp('qhat  not  positive  definite') 
end 
% 

if  abs(p'-p)>10"{-5) 
disp('p  not  synunetric') 
end 
% 

if  min  (eig(p) )<0 
disp('p  not  positive  definite') 
min(eig(p) ) 
end 
% 

%  compute  the  controller 
% 

bc=q*c ' * inv ( r2 ) ; 
cc=«-inv(  r2  )  *b' *p; 
ac=a+b*cc-bc*c+g2*q*rlinf; 
dc  =  0 ; 

% 

%compute  cl  system 
% 

[ a , b, c , d ]=lf tf (pfwd, [ 8 , 2 , 1 , 4 , 1 ] , ac,bc, cc , dc ) ; 
% 

%compute  norms 
% 

h=mynorm(a,b,c(l:2, : ) , zeros (2, 2) ) 
hinf=h(2) 

h=mynorra( a,b,c(3;4,:),zeros(2,2)) 
h2=h( 1) 
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BLACK. M  algorithum  for  solving  the  Bernstein  and  Haddad 
problem  using  the  K.  Black  technique: 

g2=l/( gam*gam) ; 
a=am; 
b»bm; 
c  =  cm; 

vl=di*dl ' ; 
v2=d2*d2 ' ; 
rl=el'*el; 
r 2  =  e2 ' *e2  ; 
rlinf=elinf'*elinf; 
sig=b*inv( r2 ) *b' ; 
sigbar=c ' *inv( v2 ) *c; 

% 

%  find  gam  =  infinity  p  and  ghat 
% 

f=a; 

g=sig; 

h=rl; 

pzero=are( f ,g,h) ; 
z=a-sig*p; 
x=sigbar ; 
qhatz=lyap( 2,x) ; 

% 

%  solve  for  q 
f*a'  ; 

g=-{g2*rlinf-sigbar) ; 
h=vl ; 

q=are( f ,g,h) ; 

% 

%  iterate  between  qhat  and  p  equations 
% 

p=pzero ; 
qhatk=qhatz ; 
er ror=l ; 
k  =  0; 

while  error>  . 000001 
% 

k=k+l 

% 

f = ( a+g2* ( q+qhatk ) *rlinf ) ; 

g=sig; 

h=r  1  ; 

pold=p ; 

p=are( f ,g,h)  ; 

% 

z  =  a-sig*p+g2*( q+qha tk/2 ) *  rlinf ; 

x=q*sigbar*q; 

qhatk=lyap( z ,x ) ; 
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% 

er ror=norm( p-pold, 'fro' ) 

% 

end 

disp( ' converged' ) 

% 

%  check  answers 
% 

if  abs(q'-q)>l0"(-5) 
disp('q  not  symetric') 
end 
% 

if  min  (eig(q) )  <  0 

disp('q  not  positive  definite') 
min(eig(q) ) 
end 
% 

if  abs(qhatk'-qhatk)>10"(-5) 
disp('qhat  not  positve  definite') 
end 
% 

if  min(eig(qhatk) )<0 
disp('qhat  not  positive  definite') 
end 
% 

if  abs(p'-p)>10"(-5) 
disp('p  not  symmetric') 
end 
% 

if  min  (eig(p) )<0 
disp('p  not  positive  definite') 
min(eig(p)) 
end 
% 

%  compute  the  controller 
% 

bc=q*c' *inv( r2)  ; 
cc*-inv( r2 ) *b' *p; 
ac=a+b*cc-bc*c+g2*q*rlinf ; 
dc=0 ; 

% 

%compute  cl  system 
% 

[a,b,c,d]=lftf (pfwd, [8,2,1,4,1] ,ac,bc,cc,dc) ; 
% 

%compute  norms 
% 

h=mynorm(a,b,c(l;2, :),zeros(2,2)) 
hinf=h(2) 

h=mynorm( a,b,c( 3 ; 4 , : ) , zeros ( 2 , 2 ) ) 
h2=h(  1 ) 


A-4 


BHD.M  algor ithum  for  solving  the  dual  of  the  Bernstein 
and  Haddad  problem: 

g2=l/( gam*gam) ; 

a=am; 

el  =  ce  ; 

c=»cy; 

dl=bw; 

d2=dyw; 

dlinf=»bd  ; 

d2inf=dyd ; 

b=bu; 

e2=deu ; 

vlinf=dlinf *dlinf ' ; 
vl=dl *dl '  ; 
v2=d2*d2 '  ; 
rl  =  el ' *el ; 
r2=e2 ' *e2 ; 
sig=b*inv( r2  )  *b' ; 
siabar=c'*inv(v2)*c, 

% 

%  solve  for  q 
% 

f=a; 

g=-(g2*vlinf -sig) ; 
h=rl; 

q=are(  f ,g,h)  ; 

% 

%initialize  p 
% 

%p=eye { n  )  ; 
p=pchest ; 

% 

%  iterate  between  qhat  and  p  equations 
% 

er ror=l ; 

while  error>.001 
% 

f =a+g2*vl inf *q-p*sigbar ; 
g=-g2*vlinf *q; 
h=q*sig*q ; 
qhat=are( f ,g,h}; 

% 

f=(a+g2*vlinf*(q+qhat) ) ' ; 

g=sigbar ; 

h=vl  ; 

pold=p; 

p=are( f ,g,h) ; 

% 

er ror=norm( p-pold, 'fro' ) 
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% 

end 

disp( 'converged' ) 

% 

%  check  answers 
% 

if  abs ( q' -q) >10* ( -5 ) 
disp('q  not  symetric' ) 
end 
% 

if  min  (eig(q) )  <  0 

disp('q  not  positive  definite') 

min(eig(q) ) 

end 

S 

if  abs ( qhat ' -qhat ) >10*  ( -5 ) 
disp('qhat  net  positve  definite') 
end 
% 

if  min( eig( qhat) )<0 
disp('qhat  not  positive  definite') 
end 
% 

if  abs ( p ' -p ) >10* ( -5 ) 
disp('p  not  symmetric') 
end 
% 

if  min  (eig(p))<0 
disp('p  not  positive  definite') 
min( eig( p) ) 
end 
% 

%  compute  the  controller 

o 

bc=-p*c ' *inv( v2 ) ; 
cc=inv( r2 ) *b' *q; 
ac=a+bc*c-b*cc+g2*vlinf *q; 
dc  =  0  ; 

% 

%compute  cl  system 
% 

[a,b,c,d]=lftf{p2, [3,2,1,1,11 ,ac,bc,cc,dc) ; 
% 

%compute  norms 
% 

h=mynorm( a,b{l;6,l),c,d); 
hinf=h(2) 

h=mynorm( a,b(l:6,2),c,d); 
h2=h(l) 
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This  is  an  example  of  an  input  file  to  the  Fortran  based 
DFP  program.  It  contains  all  dimensions,  state  space 
matrices,  numerical  parameters,  and  an  initial  controller 
guess. 


THIS  IS  THE  I. file  FOR  THE  DIRECT  METHOD 
THE  SISO  MIX  NUMBER  1  PROBLEM  GAMMA  3 . 5  MU  .000000 

THE  DIMENSIONS  I  STATE , NU , NY , ND , NE , NW , NZ 

8  1  1  2  2  2  2 

THE  PAPANETERS  GAMMA  AND  MU  (2D11.6) 

0.730D+00  O.iOOD+00 

THE  TOLERANCES  OF:  1-D  SEARCH,  CHECKSTOP  (2D11.6) 
O.lOOD-03  O.lOOD-04 


THE  A  MATRIX  ( 8F8 
-0 . 16100E+00 

.4) 

0 . lOOOOE+Ol 

0  .  OOOOOE+OO 

0. OOOOOE  +  OO 

0 . OOOOOE+OO 

0. OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

-0 . 60040E+01 

0. OOOOOE+OO 

O.IOOOOE+Ol 

0. OOOOOE  +  OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

0 . OOOOOE+OO 

0. OOOOOE  +  OO 

-0 . 58220E+00 

0. OOOOOE+OO 

0 .OOOOOE+OO 

O.IOOOOE+Ol 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0.00003E+00 

-0.99835E+01 

0. OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

0 . lOCOOE+Ol 

0 .OOOOOE+OO 

0 . OOOOOE+OO 

0. OOOOOE+OO 

-0. 40730E+00 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0 .OOOOOE+OO 

0 .lOOOOE+01 

0. OOOOOE+OO 

0. OOOOOE+OO 

-0  .  39820E+01 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

O.lOOOOE+01 

0. OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0  .OOOOOE+OO 

0 .OOOOOE+OO 

0 . OOOOOE+OO 

O.IOOOOE+Ol 

0  .OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

C. OOOOOE+OO 

0  .OOOOOE+OO 

O.OCOOOE+OO 

0 . OOOOOE+OO 

0. OOOOOE+OO 

THE  BU  MATRIX  AS 

0  .OOOOOE+OO 

BU  TRANSPOSE 

0 .OOOOOE+OO 

0.64000E-02 

0.23500E-02 

0 .71300E-01 

0 . 10002E+01 

0 . 10450E+00 

0.99550E+00 

THE  BD  MATRIX  AS 

0  .OOOOOE+OO 

BD  TRANSPOSE 

0 . OOOOOE+OO 

0 . 64000E-02 

0.23500E-02 

0 .71300E-01 

0 . 10002E+01 

0 . 10450E+00 

0.99550E+00 

0  .OOOOOE  +  OO 

0 . OOOOOE+OO 

0. OOOOOE+OO 

0. OOOOOE+OO 

0  .OOOOOE+OO 

0 .OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 

THE  BW  MATRIX  AS 

0  .OOOOOE+OO 

BW  TRANSPOSE 

0. OOOOOE+OO 

0.64000E-02 

0.23500E-02 

0.71300E-01 

0 . 10002E+01 

0.10450E+00 

0.99550E+00 

0  .  OOOOOE+OO 

0 . OOOOOE+OO 

0 . OOOOOE+OO 

0. OOOOOE+OO 

0  .OOOOOE+OO 

0 . OOOOOE+CO 

0. OOOOOE+OO 

O.COOOOE+00 

THE  CY  MATRIX 

0  .  lOOOOE+Ol 

0  .  OOOOOE-i  00 

0 . OOOOOE+OO 

0. OOOOOE+OO 

0  .OOOOOE+OO 

0 .OOOOOE+OO 

0 .OOOOOE+OO 

0. OOOOOE+OO 
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THE  CE  MATRIX 

O.OOOOOE+00  O.OOOOOE+00  O.OOOOOE+00  O.OOOOOE+00 
0.55000E-03  O.llOOOE-01  0.13200E-02  0.18000E-01 
O.OOOOOE+00  O.OOOOOe+00  O.OOOOOE+00  O.OOOOOE+00 
O.OOOOOE+00  O.OOOOOe+00  O.oOOOOE+00  O.OOOOOE+00 

THE  CZ  MATRIX 

O.OOOOOE  +  00  0 .  OOOOOE-*-00  O.OOnOOE+00  O.OOOOOE+OO 
0.5500CE-03  O.llOOOE-01  0.13200E-02  0.18000E-01 
O.OOOOOE+OO  O.OOOOOE+00  O.OOOOOE+00  O.OOOOOE+00 
O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO  O.OOOOOE+OO 


THE  DYD  MATRIX 

O.OOOObE+OO  O.lOOOOE+01 

THE  DYW  MATRIX 

O.OOOOOE+OO  O.lOOOOE+01 

THE  DE'J  MATRIX 
0 . OOOOOE+OO 

1 .  oooooE-o-; 

THE  DZU  MATRIX 
0 . OOOOOE+OO 
0 . lOOOOE+01 

THE  AC  MATRIX  (COLUMNS  1  -  4) 

-0 .170455076970+01  0 . 1 0 0 0 0 0 OOOOOD+01 

-0. 7206227137  30+01  0 . 0000  GOO  00000+00 

-0.978740782700+01  0 . 510821882940-02 

-0.158566341510+02  0 . 187567410140-02 

-0.156408137180+02  0. 5690S 7 50 3960-01 

-0.128925876180+02  0.798318823930+00 

-0.635471300420+01  0 . 83407635573D-01 

-0.339649592800+01  0 . 794567475720+00 

THE  AC  MATRIX  (COLUMNS  5  -  8) 

0 . 199 5 30’ 58900-02  0 . 399061517800-01 

0 .2  27166  499  7  90-02  0 .4  5  4  3  3  2  999580-01 

-0.248046143060-01  -0.447038340 310+00 
0 .997931422  530+0  0  -0 .2  3  3  59  55  3  5860-01 
-0.396269207180+00  -0 .637889252850+01 
-0 . 587569542620+01  -0 .109847696560+03 
-0.604225272210+00  -0. 112835454850+02 
-0. 585924577570+01  -  0.109554727540+03 

THE  BC  MATRIX 
0.154355076970+01 
0 .120222713730  +  01 
0.919532018960+01 
0 . 586950353370+01 
0.151233592580+02 
0 .736533529020+01 
0.619326642520+01 
0 .185850483400+01 


O.OOOOOOOOOOOD+00  O.OOOOOOCOOOOD+GO 
0 . 10000 00 00  0 00+01  O.OOOOOOOOOOCO-CO 
0.690417821330-02  0.993218235110-00 
0.253512793770-02  -0.654241795070-0? 
0.7691686 040 70-01  -0.19849974463D-G1 
0.10789936014D+01  -0.278456444  01  D  +  Ci? 
0.1 127  32  28  4  890+00  -0 . 2909  287982  3D-C.1 
0.107392334550+01  -0 . 277147960420+00 


0.478873821360-02  0.653009756400-01 

0  45  199599490-02  0 . 7 4  3 4  5 3 999 31D-0 

-0.‘ 55363568700-01  -0 . 7 327 8 4 7 8 3 39D-: 0 
-0 .349777559860-02  -0.386901339290-01 
-0.906542447970+00  -0.12088672144 0+02 
-0.1 247 7 369 58 <’0+0 2  -0.1799  4  8862550+0  ? 
-0. 138491440720+01  -0.17484681044  0+0  2 
-0 .13440824  04  40+02  -0 . 1 7 9 4 6 8 5 2788D+0 3 


THE  CC  MATRIX  (COLUMNS  1  -  4) 

-0.154494333890+01  0.798  159192090+00  0.1078  7  7  78  4  580+01  -  0 . 2 7 8 4 0 07 6 386D  +  0 0 

THE  CC  MATRIX  (COLUMNS  5  -  8) 

-0.589029510700+01  -0.1101  41223110  +  03  -0 . 1 3 5 1 2 5 3 3 6 5 30+02  -0  .  1 8 0 4 2 9 1  3 91 1 D+ C 3 
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Appendix  B.  Data 


Table  B.l  CENTRAL  H,  RESULTS  (II 10) 


Y 

Jredl, 

IredL 

2 . 1426 

291 . 55 

2 . 1426 

2 . 1500 

27 .601 

2 . 1500 

2 . 1600 

18 . 146 

2 . 1599 

2 . 1700 

14 . 55 

2 . 1698 

2 . 1800 

12 . 526 

2 . 1796 

2 . 1900 

11 . 190 

2 . 1894 

2 . 2000 

10.225 

2 . 1991 

2 . 2100 

9.4877 

2.2088 

2 . 2200 

8.9018 

2.2184 

2 . 2300 

8.4223 

2.2280 

2 . 2400 

8.0211 

2.2375 

2.2500 

7 .67  92 

2 . 2470 

2.2600 

7  .  3838 

2.2564 

2.2700 

7 . 1254 

2.2658 

2 . 2800 

6 . 8971 

2.2751 

2.2900 

6.6937 

2.2844 

2 . 3C00 

6.5111 

2.2936 

2 . 32 

6 . 1961 

2 .3119 

2 . 34 

5.9336 

2 . 3299 

2 . 36 

5 .7 109 

2 . 3478 

2 . 38 

5 . 5193 

2 . 3655 

2 . 40 

5 . 3527 

2 . 3830 

2 .42 

5 .2062 

2.4003 

2 .44 

5.0764 

2.4174 

2 . 46 

4.9605 

2.4343 

2 . 50 

4.7625 

2.4675 
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Table  B.l  cont 


Y 

iTedJ^ 

IredL 

2 . 5b 

4.5629 

2 . 5080 

2.60 

4.4023 

2 . 5474 

2.65 

4.2704 

2 . 5856 

2.70 

4.1602 

2.6228 

2.75 

4.0671 

2.6589 

2 . 80 

3.9  87  4 

2.6939 

2 . 85 

3 .9186 

2.7280 

2.90 

3.8586 

2.7610 

2 .95 

3.8061 

2.7931 

3 . 00 

3.7598 

2.8243 

3 . 10 

3.6821 

2 . 8840 

3 . 20 

3 .6198 

2.9403 

3 . 30 

3.5691 

2.9935 

3 .40 

3.5274 

3.0436 

3.50 

3.4926 

3,0909 

3.70 

3.4386 

3.1779 

3.90 

3.3992 

3.2555 
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Table  B.2  MIXED  RESULTS  (1110) 


Y 

WTedU 

2.1426 

2.15 

27 .6012 

2.16 

17 .9685 

2 . 17 

14 . 3460 

2 . 18 

12 .2802 

2 . 19 

10.9172 

2.20 

9.9325 

2.21 

9 . 1785 

2.22 

8 . 57  82 

2.23 

8.0863 

2 . 24 

7 .6743 

2 .25 

7.3231 

2 . 26 

7 . 0209 

2 . 27 

6 .7  57  6 

2 . 28 

6 .5193 

2.29 

6.3104 

2.30 

6.1423 

2 . 32 

5 . 8002 

2 . 34 

5.5319 

2 . 36 

5 .3050 

2 . 38 

5 .1100 

2.40 

4.9425 

2.42 

4.7945 

2.44 

4.6645 

2.46 

4.5489 

2 . 50 

4 . 3536 

2 . 55 

4 . 1601 

2.60 

4 . 007  5 

2.65 

3 . 8850 

2.70 

3.7852 

iredL 


2.1502 
2.1600 
2.1701 
2 . 1800 
2.1900 
2.2000 
2 .2100 
2.2200 
2.2300 
2.2400 
2.2500 
2.2600 
2.2700 
2.2800 
2.2900 
2.3000 
2.3200 
2.3400 
2 . 3600 
2.3800 
2.4000 
2.4200 
2.4400 
2.4600 
2 . 5000 
2 . 5500 
2.6000 
2.6500 
2 .7000 


B-3 


Table  B.2  cont 


Y 

IredI, 

iredL 

2 .1^ 

3.7030 

2.7500 

2.80 

3 .6345 

2 . 8000 

2.85 

3 . 5770 

2 . 8500 

2.90 

3.5280 

2.9000 

2.95 

3.4863 

2.9500 

3 . 0 

3.4506 

3 . 0000 

3 . 1 

3 . 3936 

3 . 0999 

3.2 

3 . 3520 

3 . 1998 

3 . 3 

3.3203 

3 . 3000 

3 .4 

3.2969 

3 .4000 

3 . 5 

3.2798 

3.4999 

3  .7 

3.2594 

3.6999 

3.9 

3.2502 

3.9273 
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Table  B.3  B  &  H  Results  (1120) 


IT..!. 

.2371 

1.7908' 

.2b74 

1 . 3401 

.211 S 

1 . 0675 

.  2977 

.  0936 

.  3179 

.  8011 

.  3382 

.7  327 

.  3b85 

.6  837 

.  37  87 

.  6464 

.  3988 

.6168 

.4188 

.  5926 

.4386 

.5722 

.  4b82 

.  5548 

.4776 

.  5397 

.4967 

.  5263 

.  5156 

.5145 

.5342 

.  5039 

.5525 

.4944 

.  57  06 

.4858 

.5883 

.  47  81 

.6058 

.  4710 

.6411 

.4583 

.647  8  ' 

.  4561 

.6488 

.4558 

.6526 

.4546 

.6593 

.4525 

.6631 

.  4513 

.6640 

.  4510 

.  6677 

4499 
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Table  B.3  cont 


.674b 

.4479 

.7119 

.4378 

.7684 

.4250 

.8220 

.4151 

.  8712 

.4075 

.9270 

.4004 

.97b0 

.  3954 

1 . 02b4 

.  3910 

1 . 0676 

.3880 

1 . 1099 

.3855 

1 . 1481 

.3836 

1 . 1839 

.3821 

1.2137 

.3811 

1.2400 

.  3804 

1.2631 

.  3799 

1.2830 

.  3795 

1 . 3002 

.  3792 

1.3140 

.3790 

1.3268 

.  37  89 

1 . 3376 

.  3788 

1.3466 

.  37  87 

1 . 3b42 

.  37  87 

1 . 3606 

.  37  86 

1 . 36b9 

.  37  86 

1 . 3704 

.  37  86 

1 . 3741 

.  3786 

1 . 377  2 

.  37  86 

1 . 3797 

.  37  86 

1.3819 

.  37  86 
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